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Abstract. The g-state Potts held theory with 2 < q < 4 in the low-temperature 
phase is considered in presence of a weak magnetic held h. In absence of the magnetic 
field, the theory is integrable, but not free at q > 2: its elementary excitations - 
the kinks - interact at small distances, and their interaction can be characterized 
^ ' by the factorizable scattering matrix which was found by Chim and Zamolodchikov. 

The magnetic field induces the long-range attraction between kinks causing their 
confinement into the bound-states. We calculate the masses of the two-kink bound 
states in the leading order in h —> ±0 expressing them in terms of the scattering 



. matrix of kinks at h = 0. 
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1. Introduction 

\C> ■ The kink topological excitations are quite common in the two dimensional field theories 
with Hamiltonian invariant under some discrete symmetry group Q. If such a symmetry 
is spontaneously broken in the ordered phase, the latter has a discrete set of degenerate 
vacua |0 a ), a = l,...,q. Then the kinks K a p, i.e. the domain walls separating two 
different vacua a and 0, behave like stable quantum particles which can propagate in 



the system. Adding a small interaction, which explicitly breaks the Hamiltonian Q- 
symmetry, lifts the degeneracy of ground states |0 Q ) and leads to confinement of kinks 
dividing the true and false vacua. This simple, but quite general scenario of confinement 
in two dimension originates to the work of McCoy and Wu [1] . Its particular realizations 
in different two-dimensional models have been the subject of a considerable interest in 
the recent years [2, 3, 4, 5, 6, 7, 8]. 

The simplest and most studied example of the two-dimensional model exhibiting 
confinement is the Ising Field Theory (IFT) [2], which is characterized by the Z 2 
symmetry group. At zero magnetic field, the theory describes free massive (apart from 
the critical point) neutral fermions, which represent in the low temperature phase the 
kinks interpolating between two degenerate vacua. At small magnetic field h > 0, 
the kinks ('quarks') become confined into pairs which form a tower of bound states 
('mesons') having zero topological charge. The meson masses M n (h) densely fill the 
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segment | [2m, oo) at h — > 0. Two asymptotic expansions describe M n {h) at small h in 
different regions of this segment. 

(i) Near the edge point 2m (i.e. for fixed n at h — >■ 0) one can use the low energy 
expansion in fractional powers of the magnetic field. Its leading term was obtained 
by McCoy and Wu [1], further corrections were found by Fonseca and Alexander 
Zamolodchikov [2, 9]. 

(ii) The masses of highly excited mesons (with n ^> 1, in particular for n ~ are 
described by the semiclassical expansion in the integer powers of h [2, 10, 11]. 

The both low-energy and semiclassical expansions for the meson masses in the IFT 
were obtained by means of the perturbative analysis of the Bethe-Salpeter equation, first 
derived for this model by Fonseca and Zamolodchikov [9]. Since their original derivation 
procedure substantially exploited the free-fermionic structure of the Ising field theory 
at h — 0, it was not applied to the models, in which kinks interact at short distances 
already in the deconfined phase. 

In this paper we address to the problem of extension of the Bethe-Salpeter approach 
to make it appropriate for calculation of meson masses in such two-dimensional models. 
The particular subject of our interest is the g-state Potts Field Theory (PFT) for 
2 < q < 4, which describes the scaling limit of the two-dimensional Potts model. In the 
latter, each site in the lattice has q different states ('colours'). At zero magnetic field, 
the Potts model is invariant under the group S g of permutations of q colours. At q — 2, 
the Potts model reduces to the Ising model. 

The zero-field PFT is integrable, i.e. it has infinite number of integrals of motion 
and the factorizable scattering matrix [12]. In the low temperature phase, the particle 
content of the (zero-field ) PFT contains q(q — 1) kinks K a p, a, /3 = l,...,q, which 
interact with each other at short distances. Their S'-matrix was found by Chim and 
Zamolodchikov [12]. 

The PFT in presence of a nonzero magnetic field h acting on one of the q colours 
has been studied recently by Delfino and Grinza [5]. Beyond other result relating to the 
T > T c phase, these authors performed qualitative analysis of the kink confinement in 
this model, classified two-kink and three kink-bound states, and conjectured evolution 
of their mass spectra with temperature and magnetic field. 

The main subject of our interest are the meson masses M n (h) in the g-state PFT in 
the low temperature phase in the limit of the weak magnetic field h — > 0. As in reference 
[5], the magnetic field is chosen acting on one colour only. We consider the mesons as 
the bound states of two kinks, which attract one another with a linear potential at large 
distances, and undergo scattering upon collisions. As the result, we express the leading 
terms in the both low-energy and semiclassical expansions for M n (h) in terms of the 
known zero-field kink S'-matrix. 

| Note, that only the mesons having masses M n (h) below the two masses of the lightest meson are the 
stable particles: M n (h) < M\(h) w 4m. Heavier mesons are the resonances since they can decay into 
a pair of light mesons. 
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The paper is organized as follows. In the next section, we remind the definition 
and some well known properties of the Potts model on the square lattice and its scaling 
limit. Sections 3 and 4 describe calculation of the meson masses in the leading order 
of the weak magnetic field in the low temperature phase at h — > +0 and at h — > — 0, 
respectively, for different values of q. In Section 5 the Bethe-Salpeter equation for the 
PFT is derived. In the IFT, the analogous Bethe-Salpeter equation provides the basis 
for the systematic theory of mesons [2]. Concluding remarks are presented in section 6. 



2. The Potts model and its scaling limit 

In this section we describe briefly the definition and few basic properties of the g-state 
Potts model in two dimensions. We start from the square lattice model and then pass 
to its scaling limit. 

Consider the two-dimensional square lattice Z 2 , and associate the discrete spin 
variable s(x) = 1,2,... ,q with each lattice site x G Z 2 . The model Hamiltonian is 
defined as 

H = -j, Yl S <x)Xy) - H ^2 S s(x), g - (1) 

<x,y> x 

Here the first summation is over the nearest neighbour pairs, T is the temperature, H is 
the external magnetic field applied along the g-th direction, and 5 a>a ' is the Kronecker 
symbol. At H — 0, Hamiltonian (1) is invariant under the permutation group S q , at 
H ^ the symmetry group reduces to S q -i. By means of mapping onto the random 
cluster model [13, 14], one can also define the g-state Potts model with noninteger values 
of q. 

The order parameters (cr a ) can be associated with the variables 

a a (x) = S s ( x ) a - -, a = l,...,q. 
q 

Parameters (a a ) are not independent, since 

5> a (s) = 0. (2) 

a=l 

The zero-field model undergoes ferromagnetic phase transition at the critical 
temperature 

Tc = iog(i + V5)' (3) 

This phase transition is first order for q > 4, and continuous for 2 < q < 4. The 

ferromagnetic low-temperature phase at zero field is g-times degenerated. For a review 

of many others known properties of the Potts model see [15, 16]. 

We shall consider only the Potts model with 2 < q < 4. In this case, the correlation 

length diverges at H — > 0, T — > T c . The Conformal Field Theory associated with this 

critical point is characterized by the central charge 

/ \ 6 . _ n . n(t — 1) , A , 

c(q) — 1 ; -, where = 2 sin — f. 4 

w *(* + !)' vy 2(t + l) V ; 
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The scaling limit of the model (1) is described by the action [5] 

A = A^ft - T J d 2 xe(i) - h J d 2 xa q (x) , (5) 

Here Wq FT corresponds to the Conformal Field Theory, which is associated with the 
the critical point. Fields e(x) (energy density) and cr q (x) (spin density) have the scaling 
dimensions 



2 V t) St(t + 1) 



Couplings r and h are proportional to the deviation of the temperature and magnetic 
field from their critical point values. 

2.1. Ordered phase in the PFT at h = 

The field theory (5) is integrable along the line h = in the (r, /i)-plane. In this paper 
only the low-temperature (r < 0) phase will be considered. At h — and r < 0, 
the S q symmetry is spontaneously broken: the model has q degenerate vacua |0 a ), 

a — 1,2, ... ,q, which are distinguished by the values of the order parameter 

v 

(crj)a = (0 a |<7 7 (a;)|0 a ) = — - (g5 7 , Q - 1) , (6) 

with some positive v. The symmetry group S q acts by permutations of these vacua. 
Elementary excitations are the q(q—l) kinks K a p{9), which interpolate between different 
vacua a and j3. Here 9 denotes the kink rapidity, which parametrizes its energy and 
momentum 

E = m cosh 9, p = m sinh 9, (7) 

with m ~ Irl 1 /! 2- ^ 9 '! being the kink mass. 

The two-kink scattering at r < 0, h = is described by the Faddeev-Zamolodchikov 
commutation relations 

K^K^) = Sll^K^Ks^), (8) 

with the scattering amplitudes 5^(^1,2), and 9 12 — 9\ — 9 2 . Due to the S q invariance, 
only four scattering amplitudes are independent, providing 

K^O^K^e*) = S (9 12 )J2K aS (9 2 )K Sf} (9 1 ) + S 1 (9 12 )K cel (9 2 )K lf} (9 1 ) , a + (5 (9) 
K^K^) = S 2 (9 12 )J2 K as(02)K Sa (9 1 ) + S 3 (9 12 )K aj (9 2 )K ia (9 1 ) , (10) 

The explicit expressions for the scattering amplitudes were determined in reference [12]: 

sinh sinh \(9 — iir) f \9\ 

01 ' sinh A (6 - 2 -f) sinh A (6 - f ) \ m ) ' 1 ' 

sinysinhy-i.) /AgX 
V ; sin f sinh A (9 - ) \m ) 
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m= nfgV (13) 



sin ^ sinh A (6 - f ) V i7r / ' 

SsW = ^n^'). (14) 



sin — \ 17T 



Parameter A is related to q as 

tt; 
"3 



Vg = 2siii^, (15) 



and 



n = BinhA(fl + if) 

\ i7r / sinh \(8 — in) 



(16) 



f°° dx sinh sinh f Q - §) 

J x smh 7^ cosh | 17T 

Note, that the values of the parameter A corresponding to integer q are: 





3 






A = 


4 


for q = 


2, 


A = 


: 1 


for q = 


3, 


A = 


3 
2 


for q = 


4. 



Zfmfc confinement in a weak magnetic field 

Application of a small magnetic field along the g-direction lifts degeneracy between the 
vacuum \0 q ) and vacua \0 a ) with a < q. In the first order in h, the shift AS between 
their energy densities reads as 

AS = 5E a - S£ q = h, for a = l,...,g-l. (18) 
Q 

It gives rise to the linear attractive potential between two kinks which interpolate 
between the stable and false vacua, 

V(x 1 ,x 2 ) = (x 2 - xi) AS (19) 

where x\ < x 2 are the spacial coordinates of the kinks. 

Depending on the sign of h, two regimes are distinguished [5]. 

• If h > 0, the vacuum |0 ? ) becomes the true ground state of the system, and the 
states |0 a ) with a ^ q become the false vacua. Magnetic field induces a long-range 
attraction between kinks leading to their confinement. Isolated kinks do not survive 
as asymptotic states of the theory, and the elementary excitations are the bound 
states of two and three kinks. § 

• If h < 0, the vacuum \0 q ) becomes metastable, and the true vacuum states |0 Q ) with 
a — 1, . . . , q — 1 are still degenerate in the energy. Elementary excitations are the 
kinks K a p(9) interpolating between the true vacua a, (3 ^ q. On the other hand, 

§ The three-kink bound states can exist if 3 < q < 4. The four-kink bound states, which are allowed 
for q = 4, should be unstable due to their decay into a pair of mesons [5] . 



Two-kink bound states in the Potts field theory 



6 



the kinks K aq {6i) and K q p(0 2 ) are confined into the bound states by the magnetic 
field. However, such bound states are unstable due to decay into isolated kinks 
K a p{9) (see [5] and the discussion below in Section 4). 

Note, the kinks in the fields theories with confinement are often called 'quarks', while 
their bound states play the role of 'mesons' (kink-antikink states), and 'baryons' (tree- 
kink states). 

3. Meson masses at h — >■ +0 

The meson mass M can be formally determined from the solution of the eigenvalue 
problem 

H(h) |tt(P)> = [E(P) + E vac ] \tt(P)), (20) 

P|7T(P))=P|7r(P)), (21) 

where T-L{h) is the Hamiltonian, and P is the total momentum operator corresponding 
to action (5), E vac is the vacuum energy, E(P) is the meson energy spectrum, which 
should have the relativistic form 

E{P) = (P 2 + M 2 ) 1 ' 2 . (22) 

Unfortunately, the explicit form of the PFT Hamiltonian 'H(h) is known only in the 
case q = 2, which corresponds to the Ising Field Theory. Below we describe briefly, 
how the meson masses at q = 2 can be calculated from the perturbative analysis of the 
Bethe-Salpeter equation in the coordinate representation [9, 10]. This procedure will be 
then naturally generalized to the case 2 < q < 4. 

3.1. Ising Field Theory case: q = 2 

In the IFT, kinks are fermions which are free at h — 0. For nonzero h, the Hamiltonian 
T-L{h) can be written as 

^-u(p)ala p -hj dxa 2 (x), (23) 

where oo(p) = \/ p 2 + m 2 is the spectrum of free fermions, and formfactors of the spin 
operator a 2 {x) are explicitly known [17]. At small h, one can treat the meson as a 
bound state of two quarks neglecting four-quark, six-quark, etc. contributions in its 
wave function. This two-quark approximation is asymptotically exact in the leading 
order in h — > 0. In this approximation, the meson mass M n can be calculated from the 
perturbative solution of the Bethe-Salpeter equation. In the coordinate representation, 
the latter equation written in the meson rest frame reads as 

2u(p) (j) {n \x) + &E\x\<P {n \x) + A£U<f)( n \x) = M n <p in \x), (24) 
where — oo < x < oo. 
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Here AS = 2/1(02) is the 'string tension', (02) is the spontaneous magnetization at zero 
field, \x\ is the distance between the two quarks, p = —id x , and <f)^ n \x) denotes the 
configuration-space wave function in the two-quark approximation 

^ n \x) = J°° ^e^(0|a ( p_ p) a p |7r„(P))p_o , (25) 

0H(-x) = -^ n \x). (26) 

The first and the second terms in the left-hand side of (24) correspond to the kinetic 
energy of two quarks, and to their long-range attraction, respectively. The linear integral 
operator U describes the the short-range interaction between quarks in the two-fermion 
sector, its kernel U(x,x') exponentially vanishes at large distances \x\ ^> mr x . 

Several perturbative schemes have been developed for the IFT Bethe-Salpeter 
equation. The most convenient for us is the procedure described in reference [10]. It is 
based on the observation, that at small h, the average distance between the quarks in 
the meson is large compared with the correlation length irT 1 . Therefore, in the right 
'transport region' at x ^> m~ l , one can very well approximate the solution 4> n (x) of the 
Bethe-Salpeter equation (24) by the bounded in the whole line —00 < x < 00 solution 
& n \x) of the equation 

[2u(p) -M n + ASx] ¥ n) (x) = 0. (27) 

After the Fourier transform this equation takes the form 

[2 u(p) - M n + i AS d p ] $ (n) (p) = 0, (28) 



providing 



where 



f(p) = 2 dqu{q) = m 2 
Jo 



' sinh 29 

9 + 



(30) 



2 

and p = m sinh 9. 

In the left 'transport region' x < 0, \x\ m _1 , function <fi( n \x) should approach 
to — $( n )(— x) due to (26). In the intermediate "scattering region" \x\ < m _1 , one can 
solve equation (24) perturbatively in h. Joining solutions obtained in these three regions 
leads to the condition which gives the meson mass spectrum M n (h). 

It turns out, that in the leading order in h, the meson masses can be obtained 
from the equation 

$(")(x)U =0 = 0, (31) 

since the odd continuous function (fi( n \x) can be approximated (in the leading order in 
h) by integral (29) in the whole positive half-axis x > 0. 

At x = 0, M n > 2 and h — > +0, the integral in (29) is determined by two saddle 
points ±p n , where 

p n = m sinh /?„ , (32) 
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and (3 n parametrizes the meson masses M n : 

M n = 2m cosh /3 n . (33) 

In the leading order in h this yields the first term of the semiclassical expansion [2, 10] 

sinh(2/3„) - 2f3 n = 2n (n - \) ( + 0(( 2 ), (34) 



$ (ri) (x 
mC 1/3 A 



= m 

x=0 



0(C) (35) 



4, 

where ( = AS/m 2 ~ h. The semiclassical expansion holds if ( <C 1 and n >> 1. 

If M n approaches 2 m, two saddle points ±p n merge at the origin, and integral (29) 
becomes proportional to the Airy function 

[°° d6 five 3 (M n — 2m) "1 "1 

(M - 2m) 

m ( 2 / 3 

in the limit h — > 0. Formulae (31), (35) give rise to the leading term of the low-energy 
expansion [1] 

— -2 = ^ n C 2/3 + 0(C 4/3 ), (36) 
m 

where — z n denotes the n-th zero of the Airy function, Ai(— z n ) = 0, n = 1,2,.... 
The low-energy expansion holds if n( < 1. It corresponds to the nonrelativistic 
approximation in the quark dispersion law oj(p) = m + p 2 /{2m) + ... in equation (27). 

3.2. Generic q 

Let us discuss now, how the procedure described above should be modified to be applied 
to the PFT with 2 < q < 4, or more generally, to models exhibiting weak confinement 
at small h, which are integrable but not free at h — 0. In this Subsection parameter q 
is allowed to take fractional values, which are assumed to be not too close to 3. 

First, it is natural to expect that the two-quark approximation can be safely used 
at small h, if the distance between quarks is much larger than the correlation length 
m _1 . We adopt two further assumptions. 

(i) At large distances x 2 — X\ ^> m -1 , interaction between two quarks K qa (xi) and 
K aq (x 2 ) forming a meson is completely described by the linear attractive potential 

V(x 1 ,x 2 ) = (x 2 -x 1 )A£ (37) 

with the string tension 

A£=-^jh + o(h). (38) 

(ii) To the linear order in the magnetic field, the quark dispersion law e(p; h) is the 
same as the free fermion spectrum u(p) = (p 2 + m 2 ) 1 ! 2 , 

e(p;h)=uj(p) + o(h). (39) 
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These natural assumptions summarize the experience gained from the IFT where 
they can be verified by means of the consistent perturbation theory based on the Bethe- 
Salpeter equation [2, 11]. In particular, radiative corrections to the quark string tension 
and dispersion law in the IFT are known to be of the third, and second order in h, 
respectively. Furthermore, the approach based on the above assumptions allows one to 
reproduce in a simple way the leading order of the exact asymptotical expansions for 
the meson masses in the IFT, see Section 2 of reference [2]. 

The Bethe-Salpeter equation (90) for the g-state PFT will be derived in Section 
5. The kernel of this integral equation is expressed in the matrix element (94) of the 
magnetization operator between the two-quark states with definite rapidities. Since the 
explicit form of this matrix element is not known for q > 2, a direct proof of statements 
(i), (ii) for general PFT is still impossible, and we shall take them as assumptions. 

A zero-momentum meson state 7r(0) can be characterized by the (q— l)-component 
wave function i/) a (x), 

4> a {x) = (Kgafa + x)K aq (xi) I 7r(0)), a = 1, . . . , q - 1, (40) 

which is well defined for large positive x ^> mT x . 

Action (5) and the true vacuum |0 9 ) are invariant under the group S q -\ of 
permutations of the first q — 1 colours. In what follows, we shall use the results of 
the S q -i symmetry analysis of the meson states, which has been done by Delfino and 
Grinza [5]. The meson states form (q — 1) multiplets n k . Here 7r fc are the eigenstates of 
the generator fl q -i of cyclic permutations of first q — 1 colours a — 1, . . . ,q — 1, 

tyf-lTT* = T^fc, (41) 

where 7 = exp[27ri/(g — 1)], k — 0,...,q — 2. 
This symmetry determines the a-dependence of the wave function ipk,a(x) for the state 

7Tfc(0), 

*Pk, a (x)=i- ka Mx)- (42) 
Due to (37) and (39), the meson wave function 4> k n \x) corresponding to the state 7rj: n ^(0) 
should satisfy equation (27) at x ^> m -1 , and therefore, we can write it in the form 

^ {X) "ft ex p{ i[/(P) A 7 Mnl +**} f - (43) 

where f(p) is given by (30). 

The further analysis depends on the value of n. 

3.2.1. Semiclassical case. For large n ^> 1 and fixed x > m -1 , the h — > asymptotics 
of integral (43) is determined by contributions of two saddle-points (32), providing 

(AC \ !/ 2 

_____ \ [^y^sinh^ + cm)e - imxsi n h ^ ) 

with C k ((3 n ) = exp 



i(& 



sinh 2/3„ ' 
2 



+ 



(44) 
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Here again the parametrization (33) is implied. Coefficients C^(/3„) and Ck(j3 n ) are just 
the in- and out- amplitudes, which characterize the quark "plane waves" which enter 
and leave at x ~ vrT x the scattering region < x < m _1 . 

At zero order in h, these amplitudes should be related by the scattering condition 

[5] 

C k (p n ) = {S 2 (2(3 n )[(q - l)4,o - 1] + S 3 (2(3 n )} C* k ((3 n ). 
The explicit form of the factor in the curly brackets in the right-hand side is: 

, (e) smhjA(| + 9)} sinh[A(i7r + fl)] for 
sinh [A(f - 9)] sinh[A(ivr - 6)] ' 

Am sinh [A(f + 0)1 
S 3 (9) - S 2 (9) = -e A ^ , - L 3 J for fc = i > ... >5 _ 



(g - 2)S 2 (0) + S 3 (0) 



2. 



sinh [A(f - 0)] 

Equations (44), (45), and (46) lead to the semiclassical quantization condition 



2tt ( // - ^ ) + U(2/3„) + 2«i(2/3 n ) + 2a 2 (2/3 r , 



sinh(2/3„)-2/3 n 
where 

ai(9) = arctan [tanh(A6>) cot(7rA)] , a 2 (9) = arctan [tanh(A0) cot (7r A/3)] . 

Equations (33), (48) determine the masses of Hq for n^$> 1. 

For the multiplet Tik , k — 1, q — 2, the analogous quantization condition 



sinh(2/3 n ) -2/3 n 



2tt I n- 1 - ) U(2 /„) + 2n 2 (2./„) 



C + 0(C 2 ) 



follows from equations (44), (45), and (47). 

For q — 4, the masses of 7r and 7r fe , /c = 1, ...,q — 2 are the same: 

r 



sinh(2/3 n )-2/3„ 



2tt I /; -- ) +iA(20 n ) 



C + 0(C) for fc = 0,..,g-2. 



(45) 



(46) 
(47) 



C + 0(C 2 ), (48) 



(49) 



(50) 



(51) 



3.2.2. Low-energy expansion Consider the leading order of the low-energy expansion 
h — > +0, n ~ 1. Functions </>&(#) are smooth in this case and vary on the scale much 
larger than the correlation length mT 1 . At large positive x ^> m _1 , they should satisfy 
the differential equation 

2m- M n - —dl + xAS <f> k (x) = 0. (52) 
m 

The regular at x — > oo solution is given by the Airy function 
<f> k (x) = Ai(t - t n ), 

where 

(M n - 2m) m 1/3 



t = (mA£) 1/3 x, t 



(A5) 2 /3 



(53) 
(54) 
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For noninteger A equations (46), (47) yield 

Vm{S 2 (0)[(q- l)4,o -1] + S 3 (6)} = -1, 

u — >0 

]im{S 3 (0) - S 2 {9)} = -1. 

9->-0 
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(55) 
(56) 



This implies the fermionic boundary condition for the wave function 0^(0) = for 
k = 0,...,q — 2 leading to the same low-energy mass spectrum (36) for all mesons 

3. 3. Weak coupling expansion for q = 3 

A separate consideration is needed at q — 3. Relations (13) - (17) reduce in this case to 
A = 1, (57) 



S 2 (0) = 



sinh(fl + i7r/3) A{0) 
sinh(# - in/3) 



Ss(9) = 0, 



, dx 2smh(x/3) , x6 

A(0) = / . , smh — • 

/ x sinnx i7r 



(58) 
(59) 



The scattering condition (45) takes now the form 

C k (P n ) = (-l) k S 2 (2p n )C* k (P n ). (60) 
In the low-energy case, the rapidities of quarks are small, 9 12 «C 1, and 

S 2 (9 1 - 9 2 ) w S 2 (0) = 1. (61) 
Thus, the boundary conditions for equation (52) should be bosonic for (f>o(x) 

O (O) = 0, (62) 
and fermionic for <f>i(x) 

0i(O) = 0. (63) 

Therefore, whereas the mass spectrum of 7Ti in the low-energy region is still given by 
equation (36), the masses of 7r are described now by 



M, 



^-2 = c 2/3 4 + o(C 



-4/3 N 



m 



(64) 



with (—z' n ) being the zeros of the first derivative of the Airy function: Ai'(—z' n ) = 0, 
n = l,2,.... 

It is tempting to associate this peculiar behaviour of the low-energy ttq mass 
spectrum in the q = 3 case with the presence of the 5-meson in the zero-field PFT 
[12]. At q = 3 and h — > +0, its mass Mb lies exactly at the edge point 2m of the 
7r-meson spectra. 

At large n>l, the zeros of the Airy function and its derivative behave as [18] 
"37r(4n-l)l 2/3 , [37r(4n - 3) n 2/3 



8 



(65) 
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Combining these asymptotics with (36), (64) and (33), we get 

( 2 ^) 3 ~ 27r L _ l\ C; for k = 0) ( 66 ) 



3! "" V 4 

(2/?n) 3 ^ _ /_ 1 



, .27r^n--JC, forfc=l (67) 
at n >> 1. 

In the semiclassical region, quantization conditions (48), (50) reduce at q = 3 to 
the form 



27r ( n ~ 7 J + 2 arctan ( tanh 2/3 " ) + U(2/3 n ) 



4/ V Vs 



1 \ / tanh 2/1, , 
2?r n - - + 2 arctan jJ-± + U(2/9 n ) 



C + 0(C 4/3 ) (68) 



C + 0(C 4/3 ) (69) 



sinh(2/3 n ) -2/3 n 
for k — 0, and 

sinh(2/3 n )-2/3 n _ ^_ 

for = 1. At /3 n — >■ these relations agree with (66), (67). 

4. Meson masses at h — >■ — 

At negative magnetic field orientated along the q-th direction, the kinks K aq and K q p 
interpolating between the true and the false vacua become confined, while the kinks 
K a p connecting two true vacua remain stable. Coupling two attracting kinks into 
bound states, one could construct the meson states both in the topological charged 
and topological neutral sectors. 

In the topological charged sector, the meson state n a/3 (0) in the 'transport' region 
x 2 — X\ > a/m, with some constant a >> 1, can be written as 

/OO l>CO 
dx 1 dx 2 \K aq (x l )K qt3 (x 2 ))iJ al 3(x2 -x x ), (70) 

oo Jxi+a/m 

where the meson wave function i]) a p{x) should satisfy equation (27) with the string 
tension 

AS = \h\((0 q \a q \0 q ) - (0>,|0 Q » = ^zr[\h\v, a ^ q. (71) 
However, the two kinks can become deconfined after the scattering process 

KaiiejKtffo) -)• K^^K^), 7 ^ q, (72) 

characterized by the amplitude So (#12) in (9). As the result, the topologically charged 
mesons n a/3 (P) are unstable already in the leading order in h for q > 3. 

On the other hand, the decay channel (72) is evidently closed for q = 3. For the 
state ^13(^1)^32(^2), the remaining scattering process 

^13(^1)^32(^2) = Sl(012)^13(0 2 )^32(0l), (73) 

for the state K 13 (9i) -^32(^2) is characterized by amplitude (12), which reduces at q = 3 
to the form 

S^O) = -e A{e) . (74) 
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Reproducing with minimal changes calculations described in Section 3, we obtain the 
masses of the topologically charged mesons 7r 12 for q = 3. For large n 1, they are 
described by the semiclassical quantization condition 



sinh(2/3 n ) -2(3 n 



2tt ( n - 7 ) + U(2/3 n ) 



C + 0(C 2 ), (75) 



4, 

and equation(33). For small n, n <^ the masses of 7r 12 are described by equation 
(36). Note, that ( = AS/m 2 ~ |fa| at negative fa. 

In the topological neutral sectors, mesons 7r aa at q ^ 4 can easily decay due to 
decoupling of kinks in the process 

K aq {0i)K qa {e 2 ) K^ejK^), 7^9. (76) 
In zero order in fa, the scattering amplitude 5*2(^12) of this channel is given by equation 
(13). 

The kink decoupling is hindered at q — 4, since ^2 (^12) vanishes in this case. As the 
result, commutation relation (14) for the mutual scattering of quarks K aq (9i), K qa (02) 
reduces at q = 4 to the form 

#a4(0l)#4a(02) = S 3 (9 12 ) K a4 {6 2 ) K Act {6 r) , (77) 

where a = 1,2, 3, and 

S 3 (0) = -e A ^. (78) 

Accordingly, the masses of 7i aa are described at q = 4 by relations (36) and (75) in the 
leading order in fa. Note, that the function A(0) depends also on parameter q through 
its (not indicated explicitly) dependence on A, see equation (17). It is natural to expect, 
that channel (76) opens in higher orders in fa making unstable all the mesons in the 
topologically neutral sector at fa — > — and q — 4. 



5. Bethe-Salpeter equation 

The heuristic approach applied in this paper is based on assumptions adopted in 
Subsection 3.2. We postulate the simple form (37) of the interaction between quarks at 
large distances, take the quark dispersion law in the form (39), and apply the boundary 
condition (45) originating from the scattering matrix at fa = 0. Though this procedure 
seems to be sufficient for determining the leading order of the meson masses at |fa| — > 0, 
it is desirable to have a more systematic theory suitable for verification of the adopted 
in Subsection 3.2 assumptions, and for calculation of higher order corrections to M n . In 
the IFT, an efficient technique based on the Bethe-Salpeter equation has been developed 
by Fonseca and Zamolodchikov [2]. In this Section we describe, how a similar Bethe- 
Salpeter equation for the PFT can be derived. We hope, that it will be used in the 
future for more consequent calculation of the weak-coupling expansion of the meson 
masses in the PFT. 
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The meson energy spectra E(P) are determined by the eigenvalue problem (20), 
(21), which we rewrite as 

/oo 
dxa q (x)\n(P)) = [E(P) + E vac ] \n(P)), (79) 
-oo 

P\tt(P)) = P\n(P)}. 

Here the Hamiltonian Pi corresponds to the integrable zero-field PFT in the ordered 
phase. Integrability of PFT at h — implies, in particular, that the Hamiltonian T-Lq 
(together with the momentum operator P) can be diagonalized by the multi-kink states 
\K«(p)): 

U \Kjp)) = [co( Pl ) + ... + co(p n )} \Kjp)), (80) 
P \Kjp)) = (p 1 + ...+ Pn ) \Kjp)), (81) 

where 

\Kjp)) = \K aihai { Pl )K aua2 {p 2 ) . . . K an _ uan (p n )) (82) 
\K a0tai (Qi)K aua2 (9 2 ) ■ ■ ■ K an _ 1:an (9 n )} 

" M0i)...w(0„)] 1/2 ' 

Oij+i (Xj, the kink momenta are ordered as oo > p\ > p 2 > ■ ■ ■ > p n > —oo, and 
9j = arcsinh(pj/m) are the corresponding rapidity variables. The kink states in the 
momentum and rapidity bases are normalized as 

(K P M\K a ^p'))=2n5{p-p'), 

{K p , a (9)\K a>p (e , )) = 2n6(e-e'). 

Note, that the particle sector of the PFT at 3 < q < 4, h = 0, contains also the 
topologically-neutral kink-antikink bound states B(p), see [12]. Of course, the mesons 
B(p) can appear together with kinks in the asymptotical in- and out-states at h — 0. 
We do not display the mesons B(pj) explicitly in formulas (80)-(83) just to avoid too 
cumbersome notations. 

Let us now turn to equations (79), concentrating on the case of a positive magnetic 
field, h > 0. Then, the meson vector \n(P)) being a topologically-neutral state in the 
sector q should admit the expansion 

|^)> = E E / 

n=2 a 1 ,...,a n - 1 ^q Jo °>Pl>->P">-°° ^> 
■\K q , ai (Pl)K ai ,a2(p2) ■ ..X an _ 1 , ? (p n ))(i ; r ? , an _ 1 (p n ) . ..if aa)ai (p2)^ai,9(Pl)k(-P)) ■ 

In complete analogy with the IFT, the two-quark approximation is based on the 
assumption, that at h — > +0, the first term with n = 2 dominates in the infinite 
sum over n in (83). Accordingly, in the two-quark approximation one replaces the exact 
eigenvalue problem (79) by its projection onto the two-quark subspace spanned by 
the basis iKq^p^K^fa)), with q and p x > p 2 , 

/oo 
dxV ( 2 q) a q (x)\n(P))= E(P) + E vac \n(P)), (84) 
-oo 

P\n(P)} = P\n(P)} 7 
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where \tt(P)) e H 2 9 \ and is the orthogonal projector on H^. Tildes distinguish 
solutions of (84) from those of the exact eigenvalue problem (79). The meson state in 
this approximation is characterized by the two-quark wave function ^/ a (pi,p 2 ), 

* Q (pi,P2) = (K q , a (p 2 )K a , q ( Pl )\7r(P)). (85) 

This relation defines \E , Q (p 1 ,p 2 ) in the domain oo > pi > p 2 > — oo. Continuation 
into the whole plane — oo < pi,p 2 < oo provided by the Faddeev-Zamolodchikov 
commutation relations (8), (10) yields 

* Q (P2,P1) = J2 S ^^ ~ ^)^/3(Pl,P2) (86) 

= s^e, - e 2 ) * a (pi,p 2 ) + 5 2 (0i - e 2 ) ^ ^(pi,p 2 ). 

/?A*,9 

Then, equation (84) takes the form 

l ~ h r°° r°° do' do' 

(p0 + u;(p 2 ) - £(P)J * Q (pi,p 2 ) = E vac V a ( Pl ,p 2 ) + -J dxj -j^- (87) 

•exp[ia;(pi + p 2 -pi -p 2 )] ^(^, Q (p 2 )K a , g (pi)| ( 7 (/ (0)|^^(p , 1 )K /3 , (/ (p 2 )) ^(pi,^). 

/3=1 

The matrix element in the integral kernel of this equation is the generalized formfactor 
of the magnetization operator in the momentum representation. It is well known, that 
such formfactors have the so-called kinematic singularities at coinciding in- and out- 
momenta of particles. Let us extract a part of these singularities which are contained 
in the disconnected 'direct propagation' terms ||, 

(^(^2)^(^)1^(0)1^^)^^)) = ^(p 2 ,Pib'i,P 2 ) (88) 
+4vr V,), [5^5( Pl - p'Mp, - p' 2 ) + S%{# x - 9' 2 )8{ Pl - p' 2 )5{p 2 - pi)] . 

Substitution of the second line of the above formula into the right-hand side of (87) 
yields 

/°° f 00 dp' dp' 
dX J (Ly 5{Pl ~ P '^ 5{J)2 ~~ Pa)*a(pi,P2) 

•exp[ix(pi +p 2 - pi -p' 2 )] = ^E vac + h(a q ) g J dx^j * a (p!,p 2 ). (89) 

In deriving of the left-hand side of (89), we have taken into account the symmetry (86). 
It is natural to expect, that the two divergent in the thermodynamic limit terms in the 
brackets in the right-hand side of (89) cancel each other. Really, in a finite system of 
the length L, the shift of the vacuum energy in presence of a positive magnetic field 
h is (in the two-quark approximation) E vac = —hL(a q ) q , and the integral f dx should 
produce the length of the system L . 

After cancelation of the infinite terms described above, the right-hand side of (87) 
becomes well defined in the thermodynamic limit, and we can safely perform in it 

|| A detailed analysis of disconnected terms in the integrals containing formfactors in integrable 
quantum field theories in a finite volume has been done by Pozsgay and Takacs [19]. 
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integration in x. Then, equation (87) takes the final form in the variables p = (pi— p 2 )/2, 
p' = (p[-p' 2 )/2: 

u(P/2 + p)+ u(P/2 - p) - E(P)\ <5> a (p; P) (90) 

30 Art' 9_1 

^^G^(pb';P) *,(;/; P). 

00 A7T p=l 

Here the meson wave function $ a (p; P) and the kernel G a p(p\p'; P) are defined as 

* a (pi,P2) = 2vr5(pi +p 2 - P)$ a (pi - P/2), (91) 
Ca^pb'; P) = ^(P/2 - p, P/2 + p\P/2 + p', P/2 - p'). (92) 

The transformation law for the function & a (p', P) under the reflection p — > — p can be 
read from (86) and (91). 

Equation (90) gives generalization of the IFT Bethe-Salpeter equation to the g-state 
PFT with 2 < q < 4. Though equation (90) and its derivation look very similar to those 
in the IFT (see Section 3 and equation (3.11) in reference [2]), an important difference 
should be pointed out. Derivation of the Bethe-Salpeter equation in the IFT is based 
on the free-fermionic basis, since the IFT describes noninteracting particles at h — 0. In 
contrast, the quarks in the PFT at q > 2 and h = strongly interact at small distances 
~ m~ l . The basis states (82), which diagonalize the PFT zero- field Hamiltonian "H > can 
be treated as its n-quark 'stationary scattering states'. Corresponding wave functions 
in the coordinate representation look like a superposition of plane waves only if the 
distances between quarks are large compared with the interaction radius. In particular, 
the wave function for the basis state \K q)a {pi)K a ^ q {p 2 )) reads as 

il>p(xi,x 2 ) = <5 a/3 exp[i(piXi +p 2 x 2 )} + S£jf(0i - 2 ) exp[i(p 2 a;i +Pix 2 )], (93) 

if x 2 — x\ ^> mr x . Here pi > p 2 is implied for the basis state, and x\ and x 2 denote 
the spacial coordinates of the left and right quarks, respectively, — oo < x 1 < x 2 < oo. 
However, in the interaction region < x 2 — x\ < m _1 , the motion corresponding to 
the state \K q ^ a {pi)K a ^ q {p 2 )) is much more complicated, and it "cannot be treated in 
terms of the wave function of a finite number of variables (because the virtual pair 
creation is possible)" [20]. Therefore, one should not understand the term 'the two- 
quark approximation' in a literal sense in the g-state PFT with q ^ 2. 

The kernel G a p{p\p'', P) of the Bethe-Salpeter equation (90) is simply related with 
the matrix element of the magnetization operator between the two-kink states in the 
rapidity basis, 

(^(02)^a lff (ei)k ff (O)|^(fli)^(^ 2 )>. (94) 

= [u(p 1 )uj{p 2 )u{p' l )uj{p' 2 )} 1/2 (X 9iQ (p 2 )X Q)9 (p 1 )|^(0)|^(p' 1 )^(p 2 )). 

The latter matrix element can be expressed with the help of the crossing relations in 
terms of the four-kink elementary formfactor of the magnetization operator 

^V««s«4(*i>fc>*sA) = ^ (95) 
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with a 4 = a . Formfactors are the central objects in the formfactor bootstrap approach 
in the two-dimensional integrable quantum field theories [17, 21, 22]. The n-particle 
formfactors are subject to a set of equations (axioms), which often allow to calculate 
them exactly (for a review see [22]). Unfortunately, the explicit expressions for all n-kink 
formfactors in the g-state PFT is known only in the Ising case q = 2. Delfino and Cardy 
[23] obtained the two-kink formfactors (0 a \o- y (0)\K a ^(6i)K l3ta (92)} of the magnetization 
operators for q — 3, 4. 

Calculation of the four-kink formfactors (95) would be crucial for determining the 
kernel G a p(p\p'; P) of equation (90). This kernel should be singular at p — ±p' due to 
the kinematic singularities of the function Gap(p2,Pi\p'i,P2)- We expect, that the leading 
singularity of G a p(p\p'; P) at p — > p' has the form 

1 1 



<W ((<7q)q ~ (<Tq) a ) 



(96) 



(p — p' + iO) 2 (p — p' — iO) 2 
This term should produce after the Fourier transform the long range attractive potential 
A£|a;| between the quarks. On the other hand, the regular at p — ±p' part of 
G a p(p\p'] P) should describe the change in their short-range interaction induced by the 
magnetic field. 

As we know from the perturbative solution of the Bethe-Salpeter equation in the 
IFT [2, 11], the regular (short-range) part of the integral kernel contributes into the 
meson masses M n {h) only in the second order in h. Physically, the additional factor h 
reflects that the two quarks bound in the meson spend at h — > almost all the time 
at large distances. They only rarely come up in the scattering region rr 2 — %i ~ in 1 
where the /i-order correction to the short-range interaction due to the regular part of 
G a p{p\p'', P) should be taken into account. 



6. Conclusion 



We extended the heuristic perturbative approach, which was originally developed 
[2, 9, 10] for calculation of the meson masses in the weak confinement regime in the 
Ising field theory, to the g-state Potts field theory with 2 < q < 4 in presence of a weak 
magnetic field h. Though the latter model is integrable at h — 0, the kinks ('quarks') 
remain to be interacting particles at h — 0. We have calculated the masses M n (h) of the 
mesons in the PFT at T < T c in the leading order in the weak magnetic field \h\ — > 
both in the low energy, and semiclassical cases. The mesons with nonzero topological 
charge were predicted for the 3-state PFT in the ordered phase at h < 0. The Bethe- 
Salpeter equation is derived for the g-state PFT with 2 < q < 4, which generalizes the 
analogous equation known in the IFT. This equation could provide a more firm basis for 
the theory, if the explicit expressions for the four-kink formfactors of the magnetization 
operator would be found. 

After the first version of this articles had appeared as a preprint, Lepori, Toth 
and Delfino [24] presented the results of their numerical investigations of the particle 
spectra in the 3-state PFT in a wide range of the magnetic fields h and temperatures 
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t by means of the Truncated Conformal Space Approach (TCSA) [25] . They confirmed 
the qualitative picture of confinement developed in reference [5], and were able to partly 
confirm our analytical predictions (36) and (64) for the low-energy part of the meson 
spectra at h — > +0. Reported in reference [24] magnetic field dependence of the masses 
of five lightest even {i = 0) and odd {i = 1) mesons at moderately large magnetic fields 
were described by formula 

M® = 2m + c®h a , (97) 

with 

c (1) 

a » 0.7, -j- » 2. (98) 
c i 

These values are in reasonable agreement with the numbers 

2 r (1) zi 

° = 3- ^f* 2 - 3 - (99) 

following from (36) and (64). However, for a complete numerical verification one needs 
to increase the accuracy of the TCSA calculations at small \h\. 

The obtained results could be developed in several directions. First, it is 
straightforward to extend them to the wide class of models exhibiting confinement, which 
are integrable but not free at zero 'magnetic field', see reference [8]. Regarding the Potts 
field theory, one can try to find from the Bethe-Salpeter equation (90) corrections to 
the meson masses at small \h\, and to study the decay mechanisms for unstable mesons 
in the higher orders in \h\. 

It is remarkable [2], that in the IFT, the Bethe-Salpeter equation reproduces with 
reasonable accuracy the mesons masses not only in the limit h — > 0, but also at finite, 
and even at large values of the magnetic field h. If this situations holds also for the 
PFT, equation (90) could be useful for nonperturbative calculations of the meson spectra 
in the PFT at finite magnetic fields. To achieve progress in all these directions, it is 
essential to find explicit expressions for the n-kink formfactors of the magnetization 
operator in the PFT. 

One more interesting open problem is to determine in the PFT the masses of 
'baryons' consisting of three quarks [5]. 

We close with the following remark. For the sake of simplicity, we calculated in 
Sections 3 and 4 the energy of a meson which has zero momentum, i.e. analyzed the 
problem in the meson rest frame. It is straightforward to modify calculations to the case 
of a generic frame, and to check (in the leading order in h), that the meson dispersion 
law E(P) really has the relativistic form (22), as one should expect. 
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